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ABSTRACT 

Using the AAOmega instrument of the Anglo- Australian Telescope, we have obtained medium- 
resolution near-infrared spectra of 10,500 stars in two-degree fields centered on the galactic globular 
clusters 47 Tuc, NGC 288, M12, M30 and M55. Radial velocities and equivalent widths of the infrared 
Ca II triplet lines have been determined to constrain cluster membership, which in turn has been used 
to study the angular extent of the clusters. From the analysis of 140-1000 member stars in each cluster, 
we do not find extended structures that go beyond the tidal radii. For three cluster we estimate a 1% 
upper limit of extra-tidal red giant branch stars. We detect systemic rotation in 47 Tuc and M55. 
Subject headings: globular clusters - Galaxy: kinematics and dynamics - Galaxy: structure - Galaxy: 
halo 
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1. INTRODUCTION 

The present structure of our Galaxy is result of a com- 
plex dynamical evolution in space and time. Globular 
clusters (GCs) are excellent tracers of this evolution be- 
cause they interact with the gravitational field of the 
Milky Way, which affects their internal dynamics (Mey- 
lan & Heggie 1997). When passing through the galactic 
disk or the bulge, tidal shocks occur and these can have 
dramatic effects on the cluster evolution (Gnedin et al. 
1999, Dehnen et al. 2004). N-body simulations have 
shown that the resulting tidal tails may be used to trace 
the orbital paths of GCs, which in turn reveals informa- 
tion about the enclosed mass of the Galaxy (e.g., Combes 
et al. 1999; Capuzzo Dolcetta et al. 2005). 

Recent star count surveys have revealed tidal tails 
stretching many degrees beyond the cluster tidal radius, 
for which the best cases include Palomar 5 and NGC 5466 
(Odenkirchen et al. 2001, 2003; Belokurov et al. 2006; 
Grillmair & Johnson 2006). These tails are believed to be 
cluster stars that have escaped due to tidal shocking. A 
crucial ingredient in these discoveries was the statistics 
of stars: in order to detect a non-circular over-density 
around a cluster, one has to eliminate galactic field stars 
very carefully. Escaped cluster members are usually iden- 
tified by their colours and apparent magnitudes, using 
optimal contrast filtering in the color-magnitude plane. 

This approach, however effective in distant and unred- 
dened clusters, cannot easily be used for nearby clus- 
ters that have strong galactic star contamination and/or 
heavy interstellar reddening. However, these clusters 
could play a very important role in constraining mod- 
els of tidal disruption, because they represent the state 
of a GC immediately before/after a gravitational shock 
by passing through the disk. The recently commissioned 
AAOmega fiber-fed multi-object spectrograph (Sharp et 
al. 2006) on the 3.9m Anglo- Australian Telescope offers 
a unique possibility for measuring radial velocities with 
1-2 km/s precision for up to 350-360 stars per expo- 
sure in two-degree field of view, which enables a spectro- 
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scopic distinction between field stars and GC members 
over several tidal radii. In this Letter we report on the 
first results of a such survey for five southern GCs. 

2. OBSERVATIONS 

The targets (Table [T]) have been selected from the cat- 
alogue of Harris (1996), using the following criteria: (i) 
nearby, thus with low or moderate reddening; (ii) small 
tidal radius (r t <20'), so that the two-degree field of the 
instrument covers an area that is several times r t across; 
(Hi) large radial velocity compared to the galactic field, 
which allows membership identification. With telescope 
time granted in August 2006, the two best GCs were M55 
and M30, both with very large radial velocities. We also 
consulted the tidal tail study of Leon et al. (2000), which 
resulted in adding NGC 288 to the targets. M12 was in- 
cluded because of the suggested severe tidal stripping (de 
Marchi et al. 2006), while 47 Tuc for its known systemic 
rotation (Meylan & Mayor 1986). Two of these objects, 
M55 and M12, lie outside the "surviving" boundaries of 
GCs in figs. 21-24 of Gnedin & Ostriker (1997), favouring 
the possibility of significant tidal tails. 

The observations were carried out on six and a half 
clear nights, between August 12 and August 18, 2006, 
with typical seeing between 1.2" and 2.2". The data 
were taken with the D1700 grating that has been opti- 
mized for recording the Ca II infrared triplet region. The 
spectra range between 8350 A and 8790 A, with a resolv- 
ing power of A/AA =10,500. Each cluster was observed 
with several field configurations, which were all centered 
on the cluster. For every configuration we exposed 60 to 
90 minutes, giving signal-to-noise ratios 50 to 250. 

Of the 392 fibers of the instrument, we used 30 for sky 
background measurements and 5 to 8 for guiding. Lists 
of target stars were prepared using the 2MASS point 
source catalog (Skrutskie et al. 2006). For M55, we se- 
lected stars that matched the color and magnitude range 
of the red giant branch. For M30 and NGC 288, we 
did not restrict the sample, while for 47 Tuc and M12 
we selected the faint end of the red giant branch. The 
K magnitudes of the selected stars fell between 8 mag 
and 15 mag. To minimize cross-talk of the fibers due 
to scattered light, each field configuration was limited to 
stars in a 3-mag wide range. For this, V- and /-band 
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Fig. 1. — Upper panels: Radial velocities and total equivalent widths of the Ca II triplet for the observed stars. Lower panels: Radial 
velocity histograms that were used to determine the velocity ranges of the clusters. Metallicity decreases from left to right. 



TABLE 1 

The observed globular clusters and their parameters. 



Cluster 


d 


[Fe/H] 


Vt 


nidal 


c 




(kpc) 




(km/s) 


(') 




NGC 104 (47 Tuc) 


4.5 


-0.76 


-18.7 


42.9 


2.03 


NGC 288 


8.8 


-1.24 


-46.6 


12.9 


0.96 


NGC 6218 (M12) 


4.9 


-1.48 


-42.2 


17.6 


1.39 


NGC 6809 (M55) 


5.3 


-1.81 


174.8 


16.3 


0.76 


NGC 7099 (M30) 


8.0 


-2.12 


-181.9 


18.3 


2.50c 



magnitudes were estimated from the 2MASS JHK mag- 
nitudes, using a set of linear transformations B,V, R, I = 
ax J + bxH + cxKs + dx, which was calibrated us- 
ing equatorial standard stars from Landolt (1983, 1992), 
cross-identified with the 2MASS catalogue (X refers to 
the Johnson bands). These transformations gave a useful 
and robust relation between Landolt and 2MASS, with 
the r.m.s. of the fits about 0.1 mag or less. 

The spectra were reduced using the standard 2dF data 
reduction pipeline that extracts the wavelength cali- 
brated spectra automatically. We performed continuum 
normalization separately and cleaned the strongest sky- 
lines that had residuals left, using linear interpolation of 
the surrounding continuum. 

In total, we collected spectra for 10,500 stars, from 
which we determined two parameters for this paper. One 
is the equivalent width of the Ca II triplet lines at 8498, 
8542 and 8662 A. Following Cole et al. (2004), we fitted 
a sum of a Gaussian and a Lorentzian function to each 
line, which was then integrated to calculate the equiva- 
lent widths and their sum EEW = Ws498+W / 8542+W / 8662- 
This parameter is sensitive to surface gravity and metal- 
licity, which can help constrain cluster membership. 
The fit also gave an initial radial velocity (RV) value, 
which was improved as follows. Using the extensive syn- 
thetic spectrum library of Munari et al. (2005), we se- 
lected best-fit synthetic spectra, degraded in resolution 
to match the observations, with \ 2 fitting, treating RV as 
a fixed parameter given by the line-profile fitting. Each 
observed spectrum was then cross-correlated with the 
best-fit model, resulting in an improved RV that was usu- 



ally within ±5 km/s of the initial value and is believed 
to be accurate within ±1-2 km/s. 

3. MEMBERSHIP DETERMINATION 

We plot EEW vs. RV in the upper panels of Fig. [TJ In 
each panel the narrow clump corresponds to cluster mem- 
bers and the broader distribution comprises field stars. 
The clusters are ordered by decreasing metallicity from 
left to right and it is also noticeable how the mean equiv- 
alent width decreases with decreasing metallicity. 

To measure the velocity ranges of the clusters, we anal- 
ysed the velocity histograms shown in Fig. [TJ For M55 
and M30, the narrow peaks of the histograms give very 
well-defined velocity limits. For the other three clusters 
we determined the velocity positions of the sudden rises 
in the histogram and applied an edge-detection algorithm 
that is searching for the extrema of the first derivative. 
The results are listed in the 2nd and 3rd columns of Table 
[2] (with uncertainties of about 1 km/s) and their straight 
averages agree within 1-2 km/s with the catalogued GC 
RVs (Table [TJ). We also experimented with fitting a sum 
of a Gaussian and a low-order polynomial background 
to measure the velocity centroid and dispersion of the 
clusters but found that Gaussians provided poor fits. 

Using the determined velocity limits, we defined the 
core samples of candidate members, i.e. stars with 
T/min <- y r <- ^m M> However, it is obvious from the 
upper panels in Fig.[TJthat even for M55 and M30, there 
is contamination from the field. We estimate the extent 
of this contamination with the following simple analysis. 

In Table [2J N tot . represents the total number of ob- 
served stars in each field, while the core samples contain 
Ny T objects. The question is how many of them have only 
by chance the same velocity as the cluster. To answer 
this, we removed the cluster peak from the histogram and 
fitted a smooth, low-order polynomial over V^. cluster ±50 
km/s. Summing up the fitted function through the clus- 
ter peak gave us an estimated number iVhist. of contam- 
inating field stars with an uncertainty of about 6-8%. 
We also counted the number of stars in the core samples 
that are located outside the tidal radius (iV rt in Table 
[2]). This number was multiplied by (60 2 /(60 2 — r 2 )) to 
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Fig. 2. — Positions of stars within the velocity range of M12, 
color coded by T,EW. The dotted circle shows the tidal radius. 



TABLE 2 

Statistics of stars in the 2-degree fields of view, 
text for the abbreviations. 



See 



Cluster 


ymm 


L^max 


AW 


AV r 


AW 


N rt 


AT' 


47 Tuc 


-32 





1696 


1149 


171 


70 


140 


NGC 288 


-51 


-38 


1219 


179 


30 


33 


35 


M12 


-52 


-31 


1780 


407 


200 


175 


192 


M55 


160 


182 


3571 


501 


8 


7 


8 


M30 


-195 


-165 


2240 


154 


7 


4 


5 



take the missing area inside the tidal radius into account 
(60' is the field radius). The resulting number, , gives 
another estimate of the field stars in the core sample. 

The good agreement between A^ist. and N' rt in most 
cases suggests that stars outside the tidal radii belong to 
the galactic field only. The largest difference was found 
for 47 Tuc, but its tidal radius covers over 2/3 of the 
whole field, so that the few stars outside do not provide 
an accurate estimate of the total contamination. 

4. CELESTIAL DISTRIBUTION OF MEMBER STARS 



Figs. [2][3] show the positions of the stars in the core 
samples, color coded by the individual EEW values (Fig. 
E} and RVs (Fig. EJ). For 47 Tuc, NGC 288 and M12 we 
show the full two-degree fields of view, while M55 and 
M30 are zoomed into the central one-degree field. The 
only evidence for non-isotropic distribution is found for 
NGC 288, where we see a NE-SW elongation that agrees 
with the reported tidal tails in the literature (Grillmair 
et al. 1995; Leon et al. 2000). However, the RVs outside 
the tidal radius do not show a coherent gradient in any 
direction, although stars in leading and trailing tails are 
expected to have systematic velocity shifts with opposite 
signs of a few km/s (Dehnen et al. 2004). On the other 
hand, there is a slight excess of stars outside the tidal 
radius, represented by the difference N' r — iVhist. = 5. 
Being positive only for NGC 288, this excess may support 
the existence of the reported tidal tails. To find a more 
conclusive result, we will use spectroscopic metallicities 
to secure membership identification for extra-tidal stars 
in a future study. For M12, the distribution of EEW 
across the field marks the cluster clearly confined to the 
tidal radius. 



Lauchner et al. (2006) discussed the relation between 
core concentration and tidal tails, showing that clusters 
with lower concentration parameter c are more likely to 
have significant extra-tidal debris, and that the cluster 
concentration is a good predictor of the presence of sig- 
nificant tidal tails. In our sample, M55 and NGC 288 
have very similarly low concentration than clusters with 
detected tidal tails (Pal 5, NGC 5466, NGC 5053), and 
yet neither of them has produced a positive detection. 

To characterize the limits of our non-detections, we 
compared the cumulative if-magnitude distributions of 
the observed stars with the whole two-degree 2MASS 
fields for M55, M30 and NGC 288. For M30, we ob- 
served 87% of all stars in the field brighter than 13.40 
mag, which is 4.7 mag fainter than the tip of the AGB. 
For NGC 288, 82% of stars brighter than 13.55 mag were 
observed (the tip of the AGB plus 5.0 mag). This means 
we can safely exclude the presence of tidal debris that 
contains 1% or more of the red giant branch stars in 
the clusters. For M55, the numbers show that 50% of 
all stars brighter than 12.9 mag that match the color- 
magnitude diagram of the cluster were observed. With 
almost 500 members captured, 1% debris that contains 
red giant branch stars would have produced two or three 
clear detections. This can be firmly excluded consider- 
ing the large SEW values of the seven extra-tidal stars 
around M55 (5.5-7.7A). We therefore conclude that the 
concentration parameter has a much weaker predictive 
power than suggested by Lauchner et al. (2006). 

To further search for any putative tidal debris, a 
matched-filter analysis was undertaken for each of the 
GC systems (Rockosi et al. 2002). Stars were selected in 
the range 17.5 > J > 8.5 and -0.5 <{J-K)< 2.2 and, 
for each GC, a background colour-magnitude distribu- 
tion was constructed by considering stars in an annulus of 
35 < r < 58 arcmins. The colour-magnitude distribution 
was also constructed using stars within the central 3.5 ar- 
cmins. The distributions were binned and smoothed to 
make them continuous, and subtracting a scaled outer 
distribution from the inner distribution clearly displayed 
the colour-magnitude sequence of the cluster itself. The 
cluster and background distributions were normalized 
and the matched-filter analysis was performed on a grid 
of 128 x 128 in a box of 40 arcmins on a side, centred 
on each cluster. Each point of the grid represented the 
centre of a box of 1 x 1 arcmin for an initial analysis, and 
5x5 arcmin for a secondary analysis, with the matched- 
filter analysis revealing what number of the stars in that 
region probably belong to the cluster. However, while 
clearly revealing the main body of each cluster, the anal- 
ysis uncovered no significant substructure that could be 
interpreted as coherent tidal debris. 

5. FUTURE WORK 

Our spectra form an excellent database for studying 
the structure and evolution of the clusters. Fig. [3] shows 
several important features that will help us to model 
their internal dynamics. Firstly, two clusters, 47 Tuc 
and M55, are characterized by large-scale asymmetry in 
velocity distribution that can be interpreted in terms of 
systemic rotation. A preliminary analysis of stellar ve- 
locities as function of position angle inferred a rotational 
velocity of 6±1 km/s and 3±1 km/s for 47 Tuc and 
M55, respectively, of which the 47 Tuc result is in ex- 
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Fig. 3. — Stellar positions for the other four clusters, color coded by the radial velocity. Note the NW-SW and E-W asymmetry in 47 Tuc 
and M55, respectively. 



cellent agreement with the 6.5 km/s obtained by Meylan 
& Mayor (1986). The influence of rotation on dynamical 
evolution may be significant, thus measuring it gives im- 
portant clues on realistic models (Combes et al. 1999). 
Our data are the first to show clear signs of rotation for 
M55. Secondly, all five clusters are characterized by large 
velocity dispersions in the center (blue and red dots sit- 
ting next to each other) that gradually goes down at the 
edge (green belts of stars halfway to the tidal radii in Fig. 
[3]). Thirdly, we do not see any evidence for increasing 
dispersion at large radii that would indicate strong tidal 
heating by the external gravitational field (Drukier et al. 
1998), the presence of surrounding dark halo (Carraro & 
Lia 2000) or a breakdown of the Newtonian dynamics in 



the weak acceleration regime (Scarpa et al. 2007). 
plan to test these theories in a future investigation. 
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